Anisotropic Scattering Rates in the t — t' — U Hubbard Model 
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We have investigated the evolution of the electronic properties of the t-t'-U Hubbard model with hole 
doping and temperature. Due to the shape of the Fermi surface, scattering from short wavelength 
spin fluctuations leads to strongly anisotropic quasi-particle scattering rates at low temperatures 
near half-filling. As a consequence, significant variations with momenta near the Fermi surface 
emerge for the spectral functions and the corresponding ARPES signals. This behavior is quite 
in contrast to the intermediate doping regime and we discuss the possible relevance of our results 
for the interpretation of photoemission spectra in cuprate superconductors at different hole doping 
levels. 
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A key issue in the efforts to understand the microscopic 
physics of high-Tc superconductors is the evolution of the 
electronic properties with doping. In recent years theo- 
retical work has continuously benefited from angular re- 
solved photoemission spectroscopy (ARPES) data for the 
electronic spectrum and the Fermi surface in the normal 
state as well as for the anisotropic energy gap in the su- 
perconducting state Q. In particular, remarkable recent 
ARPES results for the underdoped cuprates have shown 
an anisotropic normal-state pseudogap which forms be- 
low ~ 150K for the weakly underdoped materials increas- 
ing up to 300K for the heavily underdoped compounds 
with a Tc close to zero Contrary to optimally 

doped and overdoped samples the quasi-particle peak in 
the underdoped spectra is found to be very weak near the 
(tt, 0) point of the Brillouin zone (BZ) and no Fermi sur- 
face crossing is observed on the BZ boundary along the 
(tt, 0) to (tt, tt) direction. Furthermore the spin suscep- 
tibility c-axis optical |Q as well as in-plane infrared 
conductivity 1^, NMR relaxation rates j^], and inelas- 
tic neutron scattering data |lO| indicate a pseudogap in 
the low-energy excitation spectrum of underdoped com- 
pounds. Different scenarios like pair formation well above 
Tc 11 , ij] , spin-charge separation ||l3|,0 or precursor ef- 
fects near the antiferromagnetic (AF) instability [ p^p6| 
have been proposed as possible origins of these pseudogap 
phenomena. 

In this paper, we explore the combined effects of strong 
spin fluctuation scattering and Fermi surface topology 
using the 2D Hubbard model on a square lattice with a 
t — t' dispersion of the one-particle kinetic energy 



Ck — — 2t(cos kx + cos ky) — At' cos k^ cos ky 



(1) 



with nearest-neighbor (t) and next-nearest neighbor (t') 
hopping amplitudes. Near half-filling we demonstrate 
that strong quasi-particle scattering rates develop with 
decreasing temperature near the so called "hot spots" on 
the Fermi surface (FS), i.e. FS points which are con- 



nected by the AF wave vectors Q — (±7r, ±7r). As a con- 
sequence of the emerging highly anisotropic scattering 
rates the quasi-particle peaks in the spectral functions 
near the (tt, 0) points of the BZ are suppressed in com- 
parison to the momenta near hp along the BZ diagonal 
leading to highly anisotropic ARPES signals. 

Our method of choice to evaluate the renormalized 
one-particle excitations is the self-consistent and con- 
serving |]r^ fluctuation-exchange (FLEX) approximation 
]l8[ . In this approach the self-energy is given in terms 
of the spin- and density-fluctuation T-matrices Ts/(r, t) 
and Tpp{r,T) by 

j:{r,T) ^U^G{r,T)[xo{r,T) + Tpp{r,T) + T,fir,T)] (2) 

where xo(r, t) — —G{r,T)G{—r, — r) is the particle-hole 
bubble and r and r denote real space coordinates and 
imaginary time, respectively. U is the on-site Coulomb 
repulsion. The Fourier transformed Tsj(q, icjm) and 
Tpp{q,iu}m) are defined by 



Tpp{q,iuJm) 

Tsf{q,iuJrn 



1 £^Xo(q,i^m) 

2 1 + Uxoin, it^m) 

3 £^Xo(q;i^m) 

2 1 - Uxo{<i, iwm) 



(3) 
(4) 



where ujm — 2rm:T are the bosonic Matsubara frequen- 
cies at temperature T. In combination with Dyson's 



equation G 



Gn 



s, Eqs. m 



and 



form 



a self-consistent set of equations which we solve numeri- 
cally by iteration. One of the keys to the numerical solu- 
tion is the locality of (||) in space and time as well as the 
locality of (|^) and (^ in momentum and frequency space. 
This allows for a completely algebraic treatment of the 
FLEX iterations by repeated application of intermedi- 
ate Fast-Fourier-transforms p9| ] . The proper stability of 
the self-consistent cycle is achieved by solving the FLEX 
equations on a contour in the complex frequency plane 
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which is shifted off the real axis by a finite amount i'j 
with < 7 < nT/2 ||2l[. Analytic continuation to the 
real frequency axis does not encounter the usual problems 
of purely imaginary frequency methods fl^ . 

We have solved the FLEX equations on lattices with up 
to 128x128 sites using an equally spaced frequency mesh 
of 4096 points within an energy window of [— 30f, 30t]. 
The lower bound on the temperature accessible in our 
present calculations is T ~ 0.02i. This bound is set 
by the smallest width in frequency space of the AF 
paramagnon-peak in Im T^f^q, lu + 17) which can be re- 
solved for the chosen frequency mesh as well as the small- 
est momentum space width which can be treated without 
introducing finite size effects. Throughout the paper we 
will adopt an interaction strength U = At and t' = — 0.3t 
- a parameter set for which the Hubbard model exhibits 
long range AF order in the ground state at half-filling 
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FIG. 1. Im E(kF,<jJ = 0) along the Fermi line for vari- 
ous temperatures T with U = At, t' = — 0.3t and hole dop- 
ing S = 2.5%. The inset shows schematically the Fermi line 
(dashed line) and the "hot spots" (arrows) in one quarter of 
the BZ. 

We start with the discussion of our results by con- 
sidering the FS anisotropy of the one-particle scattering 
rate for small hole concentrations. Its momentum de- 
pendence is governed by two effects: First, the available 
recoil phase-space which is linked to the momentum- 
space width of the AF paramagnon peak in Im Tgf (q, us) 
and second the density of intermediate states. While 
the former quantity exhibits a variation with temper- 
ature and doping which discriminates only weakly be- 
tween the t~ and the t-t' model in the low doping limit, 
the latter quantity depends crucially on the t — t' band 
structure. This is a major source of difference between 
the one-particle renormalizations in the t- and t-t' Hub- 
bard models. In particular, we observe a non-trivial tem- 
perature and momentum dependence of the one-particle 



self-energy. This is shown in Fig. |i| which depicts the 
imaginary part of I](k, — 0) along the FS line, i.e. at 
k = kp. For T > OMt the AF peak in Im Ts/(q,u;) is a 
relatively broad structure, i.e. its HWFM in momentum 
and frequency space is tt/6 and ^ 0.2t, respectively. 
Therefore, in essence, E(ki?,w — 0) is modulated only 
by the density of states along the FS which is largest at 
the borders of the BZ and smallest at the point kp^ on 
the BZ diagonal, i.e. where kp^ = kpy. For decreasing 
temperatures the peak in Im Tsf{(i,uj) at the AF wave 
vector q = Q sharpens until we loose its accurate reso- 
lution at about T w 0.02t. This redistribution of weight 
shifts the maximum in Im l^Ckp^tu — 0) into the so- 
called "hot spots" on the FS which can be connected by 
the AF wave vector. In addition to this shift the low- 
temperature anisotropy of the self-energy is enhanced by 
roughly a factor of two. 




FIG. 2. Im G(k, ui) along three different paths in the BZ 
- schematically shown in the insets - for T = 0.03t, U — At, 
t! = -0.3f and 5 = 2.5% on a 64x64 lattice. The path 7 in 
a) is chosen along the BZ diagonal from k^. = (11, 11)^ to 
(17, 17)^. In b) the path a is parallel to the path in a) but 
runs from k^, — (21, to (32, 11)^ crossing the (hot spot) 
maximum of Im E(kF,Ci' = 0) at kfo, ~ (26,5)^. The path 
/3 in c) is along the BZ boundary from k/g = (32,0)^ to 
k^ = (32,ll)^. 

In Fig. H we show the consequence of the anisotropic 
scattering rates for the single-particle spectral function 
at T = O.OSi. Choosing a path in momentum space as 
shown in Fig. ||a) which crosses the FS at kp.^, i.e. the 
wave vector of the minimal scattering rate, a well de- 
fined, sharp quasi-particle peak is observed. In contrast, 
the quasi-particle feature is strongly suppressed in Fig. 
^) where the FS is traversed by passing through a hot 
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spot. In comparison to Fig. ga) the amplitude of the 
quasi-particle peak near is reduced by almost a factor 
of three and, moreover, it is minimal on the FS. Finally, 
choosing a path which cuts the FS on the BZ boundary as 
in Fig. pp), the quasi-particle structure is asymmetrically 
distributed as a function of momentum being more pro- 
nounced in the inverse photoemission sector. Only weak 
dispersion of the quasi-particle peak below the Fermi en- 
ergy along this cut signals a "flat-band" region close to 
the momentum (tt, 0). We note that ARPES spectra can 
be obtained from Fig. ||a)-c) by multiplying Im G(k, uj) 
with the Fermi function. Only near the momentum 
these spectra are found to display well defined quasi- 
particles which disperse through the Fermi energy, while 
the quasi-particle weight in the vicinity of the hot spots 
and the BZ boundary is substantially reduced. 




FIG. 3. Im E(kF,c<J = 0) for three points on the Fermi line 
(at maximal scattering rate, i.e. Rfq, as well as k^^, and 
kF7 as defined in the caption of Fig. 2) and hole dopings 
of 15% and 2.5% in the temperature range T = [0.02i; 0.2f]. 
The insets display the temperature dependence of the spin 
fiuctuation energy scale Usf- 

Fig. 1^ summarizes the temperature dependence of the 
imaginary part of the on-shell {lo — 0) self-energy at the 
FS for two doping concentrations S — 15% and 2.5% for 
Fermi momenta on the BZ boundary, the BZ diagonal, 
and the momentum of the maximal scattering rate, i.e. 



k^^, kpj, and \iFa, respectively. While for S = 15% the 
momenta \iFa and inpfS coincide at all temperatures de- 
picted they separate below a critical temperature T* for 
S — 2.5%. Therefore, the 'a-self-energy' in Fig.^a) con- 
tains both, an explicit temperature dependence as well 
as an implicit one which is linked to the movement with 
temperature of kpaiT) along the FS line. 

Two conclusions can be drawn from Fig. |[ First, 
we observe an intimate relation between the tempera- 
ture T*{S) and the spin fluctuation frequency uJsf{T) 
which we define by the frequency of the maximum in 
Im Ts/(q, cj) at q = Q. The insets show the tem- 
perature dependence of ujsf for both doping concentra- 
tions. Hot spots start to form once the temperature drops 
below the spin-fluctuation frequency scale, i.e. when 
T*{5) ~ uJsf{T). We find that this scenario is valid at all 
doping levels which we have investigated and moreover 
that T*{6) increases with decreasing doping concentra- 
tions. It is clearly tempting to relate T*{S) with the 
pseudogap formation temperature observed in ARPES 
jg-H and optical conductivity experiments [||. 

As a second consequence Fig. ^ suggests the existence 
of a very small energy scale which is manifest in the self- 
energy at low doping. At ^ = 15% Iml](k/,w = 0) 
clearly extrapolates to zero for vanishing temperatures 
at all FS momenta. This is consistent with Fermi-liquid 
theory (FLQ). However, at ^ = 2.5% a similar behav- 
ior can not be anticipated. This pertains to all mo- 
menta irrespective of the additional upturn of the scat- 
tering rate at kj^a below T*{S) which is partially due 
to the temperature dependent shift of ki?c(T). In or- 
der to recover low-temperature FLQ behavior we are 
forced to assume that Im T,{'kpa{T),Lj — 0,T) will ap- 
proach zero below a second characteristic temperature 
T^^Q{S) below T*{S). Since ^^/(T) is the only charac- 
teristic low-energy scale available it is natural to assume 
that r^^'3(^) - uJsfiT^^'^iS)). As is obvious from the 
inset of Fig. 0b) this temperature is expected to be very 
small for S = 2.5% and remains inaccessible within the 
numerical accuracy of our present computational scheme. 

Finally we show the quasi-particle dispersion for d — 
15% and 2.5% in Fig. |. The quasi-particle energies 
have been determined from the zeroes of the real part 
of the inverse Green function, i.e. from -^(k) — e(k) — 
Re S(k, £^(k)). Far off the Fermi energy this figure seems 
to suggest a nearly rigid-band picture. However, close to 
the FS the situation is less trivial, in particular close to 
the wave vector k = (tt, 0), where quasi-particle energies 
get strongly renormalized. With decreasing doping con- 
centrations we find the flat band region around this wave 
vector to be 'pinned' to the Fermi energy and - surpris- 
ingly - deformed towards a stronger dispersion. While 
for S — 15% the quasi-particle energy i?(7r, 0) ~ — 1.09t 
is very close to the bare energy of e(7r,0) = — 1.20t, 
the pinning gives rise to a substantial renormalization at 
S = 2.5% where we find E{tt, 0) ~ -0.73i. Qualitatively 
similar pinning, however quantitatively less pronounced, 
can be observed along the BZ diagonal. 
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FIG. 4. Quasi-particle dispersion relation for two different 
hole doping concentrations along the closed path in the BZ 
(0,0) ^ {tt,tv) (7r,0) (0,0). The horizontal lines indi- 
cate the chemical potentials for both doping levels. The inset 
shows a magnification of the vicinity of the (7r,0) point. 

Focusing on the low-energy sector in the vicinity of the 
Fermi surface the inset in Fig. ^ details yet another re- 
markable, though weak effect. While a continuous band 
dispersion is evident from this inset, with a well defined 
Fermi surface crossing on the BZ boundary, the quasi- 
particle energy develops a step like feature around (tt, 0) 
for 5 = 2.5%. This is a precursor of the energy gap open- 
ing at half-filling and induces a pseudogap structure in 
the density of states. We have verified this feature to 
be robust against finite size scaling on up to 128 x 128 
lattices. Clearly the magnitude of the pseudogap is too 
small to account for the recent corresponding findings of 
ARPES |-|]. A tempting remedy for this deficiency is 
to decrease the ratio of \t' /t\ and thereby enhance the 
effects of nesting. Indeed, within our FLEX scheme, we 
found that reducing t' leads to a pronounced increase of 
the pseudogap structure and moreover to the appearance 
of a split-off precursor band similar to that observed in 
ARPES H. This is also consistent with QMC data ob- 
tained on the pure t-model However, we believe that 
this route to pseudogap behavior is inconsistent, both 
with the Fermi surface topology of the cuprate materi- 
als as well as with the strong anisotropy of the quasi- 
particle scattering rate, which both require the inclusion 
of a finite hopping amplitude t'. Therefore a proper un- 
derstanding of the combined effects of anisotropy and 
pseudogap formation requires further investigations. 

In conclusion we have studied the anisotropic elec- 
tronic properties of the two-dimensional U — t — t' Hub- 
bard model which result from the unique combination 
of the shape of the Fermi surface and strong scatter- 
ing from spin fluctuations. With decreasing hole dop- 
ing and lowering of the temperature hot spots develop 
on the Fermi surface with significantly enhanced scatter- 
ing rates offering a physical explanation for the observed 
strong anisotropics of the ARPES signals in the under- 



doped cuprates. For realistic values of the ratio t' /t clear 
Fermi surface crossings are still maintained in our FLEX 
calculation and therefore a pronounced pseudogap forma- 
tion is not obtained. An extension to bilayer type Hub- 
bard models is a possible path for further study since the 
bilayer splitting into a bonding and an antibonding band 
provides another important low energy scale which may 
be of particular relevance near the (tt, 0) point of the BZ. 
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